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Abstract 

We calculate the penguin form factor for the B r]i~^i~ decay. This form factor is 
calculated in light cone QCD sum rules, including contributions from wave functions 
up to twist-4, as well as mass corrections of the light rj meson. 



PACS numbers: 13.20.He, 13.20.-v, 14.40.-n 



*e-niail: taliev@metu.edu.tr 
^e-mail: savci@metu.edu.tr 



1 Introduction 

Rare B meson decays, induced by flavor changing neutral current (FCNC) b — > s{d) transi- 
tion, provide potentiaUy the stringiest testing ground for the Standard Model (SM) at loop 
level. These decays are also very suitable looking for new physics beyond the SM. Among all 
decays of B mesons, the semileptonic decays receive special attention, since their study offer 
one of the most efficient ways in determination of the Cabibbo-Kobayashi-Maskawa (CKM) 
matrix elements. From experimental side, there scheduled an impressive program for study 
of both inclusive and exclusive B-decays in B factories, BaBar and Belle, as well as LHC-b 
machines. CLEO Collaboration [1] has measured the branching ratios of TT~i'^u and 

B p-i+u decays, from which it is obtained that |Kb| = (3.25 ± 0.14t[|J^ ± 0.55) x 10"^ 
In extraction of \Vub\ from B — > 7r(p)£z/ decay, main theoretical uncertainties come from 
B — > 7r(p) transition form factors. For an accurate calculation of the CKM matrix elements, 
hadronic form factors need to be determined more reliably. 

It should be noted that the decay modes of i? ^ K£^i~ {£ = e, /x) has recently been 
observed with B{B Ki+i") = (0.75l[jjf ± 0.09) x lO'^ [2] and (0.78tl]J^t[]:l^) x 10-^ 
[3, 4]. At BaBar, an excess of events over background with 2.8o" has been observed for the 
B K*i+i- decay with B{B K*i+i-) = {l.68°Jo%s ± 0.28) x lO"*^ [4]. 

In this work we calculate the penguin form factor of the B — > rj£'^£~ decay in light cone 
QCD sum rules. The form factors induced by the vector current in — > rjii' decay has 
already been calculated in light cone QCD sum rules in [5]. It should be mentioned here 
that B ^ 1] form factors are related to the B ^ tt form factors through SU(3) symmetry, 
which are calculated in light cone QCD sum rules in [6] . A detailed description of the light 
cone QCD sum rule and its applications can be found in [7, 8]. 

Interest to i? — >■ r]t^£~ and B — >■ vlt^l~ has its grounds in the fact that they can 
give information about 77-77' mixing angle [9, 10]. Soon B factories will provide much more 
data and therefore a more reliable determination of the transition form factors and as a 
result a more precise determination of will be possible. The extraction of from 
the B r\{r]'^t^t~ decay would present an efficient and complementary alternative to its 
determination from B — > 'k{^p)1'^1~ decay. 

The present work is organized as follows. In section 2, we calculate the sum rule for 
the penguin form factor of the B — > 'r]l'^l~ decay. Section 3 is devoted to the numerical 
analysis and the conclusion. 



2 Light cone QCD sum rules for the penguin form 
factors 'uv B ^ r\ transition 

The penguin form factor of the B^ ^ rj transition is defined as 



(1) 



TUB + rrir, 

The starting point for the calculation of the form factor Jt in Eq. (1) is the following 
correlator function: 

(l) = ij d'xe''^^ {v{p) \T{qa^,q''{l + -f,)b{x)b{0)i{l - 75)?} | 
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which is calculated in an expansion around the light cone x"^ — 0. The main reason for 
choosing the chiral — 75)5 current instead of the bi'j^q current which has been used 
in the calculation of the B ^ tt form factor [8], is because twist-3 wave functions do not 
contribute for this choice, which are the main inputs of the light cone QCD sum rules and 
which bring about the main uncertainty to the results [11]. 

Following the general idea QCD sum rules to obtain the penguin form factor is by 
matching the representation of the correlator function in hadronic and quark-gluon lan- 
guages. Let us first consider the hadronic representation of the correlator function. By 
inserting a complete set of states with the same quantum numbers of the B meson between 
the currents in the correlator, and singling out the pole term of the lowest pseudoscalar B 
meson, we get 
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where the sum in Eq. (3) describes the contributions of the higher states and continuum. 
For the invariant amplitude 11^ one can write a general dispersion relation in the B meson 
momentum squared {p-\- qY as 



U^{q\{p + qf)= jds- 
The spectral density corresponding to (3), is 



(p + qf 
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where we have used the definition 
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The first term in Eq. (5) represents the ground state B meson contribution and p^{s) 
corresponds to the spectral density of the higher resonances and the continuum. The 
spectral density p^{s) can be approximated by invoking the quark-hadron duality ansatz 



(6) 



where sq is the continuum threshold. As a result, the hadronic representation of the invari- 
ant amplitude 11^ takes the following form 



mq^Wsfs 
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-|- subtractions 
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In order to obtain the sum rule for fj-iff)-, we proceed to calculate of the correlator 
function from QCD side. This can be done by using the light cone OPE method. For this 
purpose, we work in the large space-like momentum regions {p + qY — ml <^ for the 
hq channel and (f' <^ — 0(few GeV"^) for the momentum transfer, which correspond 
to the small light cone distance x and are required by the validity of the OPE. After 
contracting h quark field, we get 



n^(P, (l) = ^j t^'^e*"" {V{P) \q{x)a^,q%l + l^)S\x, 0)z(l - 75)? 



(8) 



where S^{x, 0) is the full quark propagator. In presence of the background gluon field, its 
explicit expression can be written as 
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where the first term on the right hand side corresponds to the free quark propagator, G"^ is 
the gluonic field strength and Qs id the strong coupling constant. We see from Eqs. (8) and 
(9) that, in order to calculate the theoretical part of the correlator, the matrix elements of 
the nonlocal operators between 1] meson and vacuum states are needed. 

Here we would like to remark that in the following calculation rj-r]' mixing will be 
neglected, since in octet -singlet basis this angle is about 9 ~ 10° [12]. Hence, in the 
above-mentioned basis, the interpolating current for 77 meson is chosen as the SU(3) octet 
axial-vector current 
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In order to simplify the notation we will use qVq to denote 

= (wr^^i + dV^d - sP^s) , 

and introduce = /^/ v^. Here, is the leptonic decay constant of 77 meson and is to be 
determined from the relation 



(11) 

It is easy to see from Eqs. (8) and (9) that the terms containing even number of Dirac 
matrices do not give any contribution. Remaining matrix elements can be parametrized in 
terms of r] meson functions up to twist-4 defined as 

/■I r 1 
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where 
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In Eqs. (12)-(14), the function <^^(m) is the leading twist-2, A{u)^ (f\\{ai)^ ^Acxi), (^j|(Q;j), 
and '^xio-i) are all twist-4 wave functions. Inserting Eqs. (12)-(14) and Eq. (9) into Eq. 
(8) and completing integration over the variables x and A;, we get for the invariant structure 
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The next and the last step in obtaining the sum rule for penguin form factor is to carry out 
the Borel transformation with respect to the variable (p + g)^ which enhances the ground 
state contribution and suppresses contributions of the higher states and the continuum. 
Finally, matching this result with the corresponding invariant amplitude that is calculated 
in hadronic and quark languages, we get the sum rule. Subtraction of the continuum contri- 
bution is performed by using quark-hadron duality (more about subtraction of continuum 
and higher state contributions in hght cone QCD can be found in [14, 15]). Performing 
Borel transformation in Eq. (15), we get for the theoretical part 
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In the same manner, performing Borel transformation in Eq. (7) and equating it to Eq. 

(16) , we finally get the following sum rule for the penguin form factor 

/^(g^) = e-l/^^ (n^)"" . (17) 

3 Numerical analysis 

In this section we present the result of our numerical calculations on penguin form factor 
fxiff)- It follows from Eqs. (16) and (17) that the main input parameters of the sum rule 

(17) are the r] meson wave functions. The explicit expressions of the wave functions ^t^{u), 

A{'u), (/9[[(aj), ip^{a,j), and (^j_(aj) are all given in [13]. The other necessary input 

parameter of the sum rule is the leptonic decay constant F^. As has already been noted, 
we will r]-ri' mixing. Furthermore, since 77 meson is an isoscalar, we have 

Tpd ciu = p f ^ 

where for the leptonic decay constant rj meson, we quote the result of a recent analysis 
which predicts frj = 159 MeV [16]. Moreover, the leptonic decay constant B meson is 
chosen to have the value Jb = 160 MeV [14, 17]. 

Having all these input parameters at hand, we proceed carrying out numerical calcula- 
tions. First of all, since is an auxihary Borel parameter, we must find a region of 
where a physically measurable quantity be practically independent of it. The lower bound 
of is determined by the fact that nonperturbative terms must be subdominant. The 
upper limit of is determined by the condition that the higher states and continuum 
contributions are less than, for example, 30% of the total result. Our numerical analysis 
shows that both conditions are satisfied in the region 8 GeV'^ < < 16 GeV'^. Moreover, 
it should be emphasized that light cone QCD sum rule predictions are reliable in the region 
of momentum transfer square, i.e., < ml — 2mfeA, where A is a typical hadronic scale 
having the value A ~ 0.5 GeV, which yields ^ 18 GeV^. 

In Fig. (1) we present the dependence of the form factor fTil^) on the Borel parameter 
at different values of momentum transfer square, q^ — GeV'^, q^ — b GeV"^ and 
(f — 10 GeV^^ at two different choices of the continuum threshold Sq = 35 GeV^ and 
So = 40 GeV^. We observe from this figure that, fj, seems to be practically independent of 
the Borel parameter M^, as varies in the region 8 GeV^ < < 16 GeV^. 

Having this window for M^, we next study the dependence fT^I^) on g^, at three fixed 
values of the Borel parameter = 8 GeV^, = 12 GeV^ and = 16 GeV^, picked 
obviously from the above-mentioned working region of M^, again at two fixed values of the 
continuum threshold, sq = 35 GeV^ and sq = 40 GeV^, as before. Depicted in Fig. (2) is 
the dependence of the form factor on the momentum transfer g^, which clearly demonstrates 
that /J(0) — 0.16 ± 0.03. As we have noted earlier, the prediction by the hght cone QCD 
sum rule is not rehable in the region q^ >!'& GeV^. In order to extend the present result 
to whole physical region, we look for some convenient parametrization of the form factor 
in such a way that in the region 4m| < < 18 GeV^ this parametrization coincides with 
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the light cone QCD sum rule prediction. The best parametrization of with respect to 
can be written in terms of three parameters in the following way 



fm 



(18) 
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1 — ap — 9- + bp — 9- 



For the values of these parameters for the penguin form factor we obtain ap — 1-08 and 
bp — 0.09, where the quoted errors can be attributed to the variation in sq and M^. As 
has already mentioned earlier, the form factor for the B rj transition can be related to 
the corresponding B ^ rr transition form factor through SU(3) symmetry. For example, 
the value /^(g^ = 0) = 0.17 is obtained using SU(3) symmetry seems to be in quite a 
good agreement with our prediction of JtIq^ = 0). In conclusion, we have calculated the 
penguin form factor for the B — > rji'^i" decay in light cone QCD sum rule method, including 
contributions of wave functions up to twist-4 and mass correction of the 77 meson. 
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Figure captions 



Fig. (1) The dependence of the form factor fj. on the Borel parameter at 5^ = GeV'^, 
5 GeV'^, and 10 GeV"^, at fixed values of the momentum threshold sq = 35 GeV'^ and 
So = 40 GeV^. 

Fig. (2) The dependence of the form factor fj- on the momentum transfer at = 
8 GeV'^, 12 GeV'^, and 16 GeV"^, at fixed values of the momentum threshold sq = 35 GeV"^ 
and So = 40 GeV^. 
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